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Abstract: A variety of supersymmetric models give rise to a split mass spectrum char¬ 
acterized by very heavy scalars but sub-TeV gauginos, usually with a wino-like LSP. Such 
models predict a thermally-produced underabundance of wino-like WIMP dark matter 
so that non-thermal DM production mechanisms are necessary. We examine the case 
where theories with a wino-like LSP are augmented by a Peccei-Quinn sector including an 
axion-axino-saxion supermultiplet in either the SUSY KSVZ or SUSY DFSZ models and 
with/without saxion decays to axions/axinos. We show allowed ranges of PQ breaking 
scale f a for various cases which are generated by solving the necessary coupled Boltzmann 
equations. We also present results for a model with radiatively-driven naturalness but with 
a wino-like LSP. 
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1. Introduction 


Supersymmetric models with anomaly-mediated SUSY breaking [1] (AMSB) provided a 
strong raison d’etre for considering the case of a wino-like lightest SUSY particle, or LSP. 
Originally, such models were built with a “sequestered”- rather than a hidden- SUSY 
breaking sector. The sequestered sector could be located on a brane which was separated 
from the visible sector brane in an extra dimensional space-time. In such a case, tree level 
supergravity contributions to soft SUSY breaking terms were absent and the dominant 
contribution to soft terms came from the superconformal anomaly. Since the soft terms 
were all of order m so ft ~ m 3 / 2 /( 167 r 2 ), then values of gravitino mass m 3 / 2 ~ 30 — 100 TeV 
were required to generate a weak-scale sparticle mass spectrum. The weak-scale gaugino 
masses were expected to occur in the ratio : M 2 : M 3 ~ 3 : 1 : 8 , resulting in a wino- 
like LSP as the dark matter candidate. The thermally-produced relic density of a wino-like 
LSP is typically [2, 3] 

fl—’/i 2 ~ 0.12 (M 2 /2.5 TeV ) 2 . (1.1) 

The measured dark matter abundance Ucdm ^ 2 = 0.12 is then saturated for a wino of mass 
m Wi - 2.5 TeV. For lighter winos, non-thermal production mechanisms such as 

WIMP production from moduli decay were invoked [4]. 

While the simplest AMSB models provided solutions to the SUSY flavor, CP and 
gravitino problems, they retain the problem of predicting tachyonic slepton masses. More 
recently, they may have fallen into disfavor due to the discovery [5, 6 ] of the Higgs boson 
with mass rrih = 125.5 ± 0.5 GeV. In the minimal AMSB model, this value of Higgs 
mass requires m 3/ / 2 ~ 1000 TeV so that the sparticle mass spectrum lies in the multi- 
TeV region which seems to seriously compromise even the most conservative measures of 
naturalness [7, 8 ]. 

Even well-before the Higgs discovery, related models with a wino-like LSP were emerg¬ 
ing. These include 

• PeV SUSY [2, 9], 

• split SUSY [10, 11, 12], 

• G2MSSM [13], 

• models with strong moduli stabilization [14], 

• pure gravity mediation [15, 16] and 

• spread SUSY [17, 18]. 

These models differ from the original rnAMSB model in that they predict a split 
spectrum with scalars ranging from 25 TeV all the way to 10 8 TeV- well beyond the 
reach of collider experiments. In contrast, the gauginos typically lie in the 0.1 — 3 TeV 
region so that the lower range of values would be accessible to LHC searches. In most 
of these models, the gauginos adopt either the AMSB-form [9, 15] or a mixed anomaly 
plus loop contribution form [13, 17, 18] which also typically gives rise to a wino-like LSP. 
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The SUSY fj, parameter is variable between these several models and may be as small 
as ~ 1 TeV [13, 16] or as high as hundreds of TeV [15]. While the predicted thermal 
abundance of wino-like WIMPs saturates the measured value for a wino mass of ~ 2.5 
TeV (so the gaugino spectrum would be well beyond reach of LHC), for lower M 2 values 
a thermal underabundance of WIMPs is expected and some non-thermal DM production 
mechanism is needed. Usually, this has involved some form of moduli production and 
decay [4, 19, 20, 21] (for recent reviews, see Ref’s [22, 23]). 

In the present paper, we instead look at non-thermal wino production from the Peccei- 
Quinn (PQ) sector. 1 By invoking a PQ sector in supersymmetric models [25] the axion 
supermultiplet also contains an R-parity-even spin-0 saxion s and an R-parity-odd spin-1/2 
axino a. This approach has several advantages: 

• a PQ sector seems necessary to solve the strong CP problem in the QCD sector [26, 
27, 28], 

• invoking PQ charges for Higgs multiplets offers a means to forbid the appearance of 
a Planck scale /i term while re-generating a weak-scale term (solution to the SUSY 
fi problem) [29], 

• while the presence of the PQ sector can act to augment the wino abundance- for 
instance by axino and/or saxion decays- the axion abundance can always be adjusted 
to make up any remaining DM abundance which may be needed. 

To explore this situation, we will adopt a benchmark model which encapsulates the 
dark matter physics expected in the above list of models. This benchmark point- labelled 
as CSB for “charged SUSY breaking” [9]- contains scalar masses around the 72 TeV region 
while gauginos lie in the 0.2 — 2 TeV range. The thermally-produced WIMP abundance 
is predicted to be O^/i 2 ~ 0.002- a factor ~ 60 below the measured value. Such a 
low thermal WIMP abundance requires additional dark matter production mechanisms to 
match experiment. In the case presented here, the dark matter is actually composed of 
both WIMPs and axions. While WIMPs can be produced thermally, they can also be 
produced via axino, saxion and gravitino production and decay in the early universe. In 
addition, saxions produced via coherent oscillations (CO) can inject late-time entropy into 
the early universe, thus diluting any relics already present. Axions can be produced as 
usual via CO [30, 31], but can also be produced thermally and via saxion decay. 

While the models listed above are motivated by a variety of theoretical and phenomeno¬ 
logical considerations, we note that collectively the entire set is highly fine-tuned in the 
electroweak sector, since the weak scale values of vn 2 H and /i 2 would have to be adjusted to 
very high precision to gain a Z mass of just 91.2 GeV. Thus, for contrast, we also examine 
a SUSY model with radiatively-driven naturalness [32] but with a wino-like LSP [33] with 
fine-tuning at just the 10% level (labelled as RNSw). 

In Sec. 2, we briefly review a variety of models with split spectra and a wino-like LSP. 
We also present a SUSY model with radiatively-driven naturalness and a wino-like LSP 

1 An earlier look at non-thermal production of winos in AMSB models was given in Ref. [24]. 
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for comparison. In Sec. 3, we briefly review our coupled-Boltzmann equation evaluation 
of mixed axion/wino dark matter (more details can be found in Ref. [34]). In Sec. 4, we 
present the results of our coupled Boltzmann computation of the mixed axion/wino dark 
matter abundance in the CSB and RNSw benchmark models. In Sec. 5, we expand our two 
benchmark points to model lines to examine how our results depend on the SUSY mass 
spectrum. Our overall conclusions and a summary plot are given in Sec. 6. 


2. Survey of some models with a wino-like LSP 


2.1 PeV SUSY 


In Ref. [2, 9], it is argued that the PeV scale (with m(scalars) ~ 771 , 3/2 ~ 1 PeV=1000 TeV) 
is motivated by considerations of wino dark matter and neutrino mass while providing 
a decoupling solution [35] to the SUSY flavor, CP, proton decay and gravitino/moduli 
problems. This model invoked “charged SUSY breaking” (CSB) where the hidden sector 
superfield S is charged under some unspecified symmetry. In such a case, the scalars gain 
masses via 




4fs. 


while gaugino masses, usually obtained via gravity-mediation as 


( 2 . 1 ) 


d 2 e-^-WW => -^-AA, 


M P 


M P 


( 2 . 2 ) 


are now forbidden. Then the dominant contribution to gaugino masses comes from AMSB: 

(2.3) 


33 o 2 

Ml = TI^2 m3/2 ~ m 3/2/ 12 0, 


02 


m 2 = ^ 2 ^ 3/2 ~ rn 3/2 / 360, 


03 


M 3 = ~ 3 ][^2 m 3/ 2 ~ -"*3/2/40. 


(2.4) 

(2.5) 

( 2 . 6 ) 


Saturating the measured dark matter abundance with thermally-produced winos requires 
m w ~ M 2 rs_/ 2.5 TeV which in turn requires the gravitino and scalar masses to occur at 
the ~ 1000 TeV (or 1 PeV) level. The author remains agnostic as to the magnitude of //, 
although fi M 2 is expected. 


2.2 Split SUSY 

In Split SUSY [10, 11, 12, 36, 37], SUSY is still required for gauge coupling unification and 
for a dark matter candidate, but naturalness is eschewed in favor of a multi-verse solution 
to the gauge hierarchy problem. In such a case, matter scalars can exist with masses 
typically at some intermediate scale ~ 10 8 TeV while SUSY fermions (gauginos and 
higgsinos) are protected by chiral symmetry and can be much lighter. Split SUSY can 
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be realized under charged SUSY breaking as in PeV-SUSY or via Scherk-Schwartz SUSY 
breaking in extra dimensions [10]. Here, one might expect 

m(gauginos) ~ m(higgsinos) <C m(scalars) (2-7) 

where the authors remain agnostic concerning whether the wino or bino might be lighter. 
Typically, binos should overproduce dark matter so that a wino/higgsino admixture might 
be expected. 

2.3 G2MSSM 

In string/M-theory models which are compactified on a manifold of G 2 holonomy [13], 
one expects a gravitino mass m 3 / 2 ~ 25 — 100 TeV along with a cosmologically relevant 
moduli field with similar mass [38]. The matter scalar masses are of order ~ m 3 / 2 but 
gaugino masses can be much lighter. Typically, a wino LSP is to be preferred [20]. The 
superpotential \i term is generated with value ~ 1 TeV so that these models tend to be 
more electroweak-natural than split SUSY. 

2.4 Models with strong moduli stabilization (Kallosh-Linde or KL) 

In string theory, an outstanding problem exists in the need for vacuum stabilization of 
moduli fields. In the KKLT construction [39], one constructs a stable supersymmetric anti- 
deSitter vacuum, but then uplifts to a deSitter vacuum via SUSY breaking. In KKLT, the 
volume modulus mass m a is expected to be comparable to the gravitino mass m 3 / 2 . These 
models give rise to soft SUSY breaking terms characterized by comparable moduli- and 
anomaly-mediated contributions [40]. However, these models suffer from vacuum destabi¬ 
lization during inflation unless the Hubble constant H < m 3 / 2 . Such inflationary models, 
while possible, are often unwieldy and inelegant [14]. 

An alternative approach known as strong vacuum stabilization invokes instead a race¬ 
track superpotential for the volume modulus, leading to a far heavier modulus mass 
m a ~ 10 15 GeV and allowing for vacuum stability in models of chaotic inflation [14]. 
In this Kallosh-Linde (KL) case [41], the soft SUSY breaking scalar masses are comparable 
to m 3 / 2 , but the gaugino and trilinear soft terms are suppressed by a factor of m 3 / 2 /m a . 
The dominant contribution to gaugino masses comes from anomaly-mediation. Requiring 
a wino LSP without too much relic density then fixes m 3/ / 2 < 1000 TeV. Thus, one gains 
a model of split SUSY with PeV-scale scalar masses but with TeV-scale gauginos with an 
AMSB mass pattern. The /x parameter is also expected to be ~ m 3/ / 2 [42] so a high degree 
of electroweak fine-tuning is needed. 

2.5 Pure gravity-mediation 

In pure gravity mediation (PGM) models [15], it is assumed that matter scalar masses 
are developed at tree level and so have masses m~^ H ~ m 3 / 2 ~ 1000 TeV while gaugino 
masses are suppressed since no SUSY breaking fields are assumed to be singlets under any 
symmetries. The gaugino masses arise via anomaly mediation so the wino is expected to 
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be the LSP. The /r term and SUSY breaking bilinear B are also expected to be at the 7713/2 
scale leading to 

m g,b,W < m qlH.h ~ 100 TeV ( PGM ) ( 2 - 8 ) 

although a recent incarnation also allows for light higgsinos [16]. 

2.6 Spread SUSY 

In Spread SUSY [17, 18], additional spatial dimensions are assumed so that the 4-d reduced 
Planck scale Mp is enhanced by a volume factor over the fundamental scale M*. Then, 
if the hidden sector SUSY breaking field X is charged under some symmetry, gaugino 
masses are generated only via anomaly- mediation while scalar masses are generated via 
gravity-mediation. One expects a mildly split- or spread- SUSY spectrum characterized 
by 

m w,b,g < m 3/2 ~ rn- h < m- iH (2.9) 

where the wino is the LSP with sub-TeV masses and the matter scalar masses may lie in 
the 10 2 — 10 3 TeV range while the higgsinos are intermediate between these two. 

2.7 Natural SUSY with wino-like LSP 

In SUSY with radiatively-driven naturalness [32], the W,Z,h mass scale arises naturally 
due to a supersymmetric ^ parameter with /i ~ 100 — 300 GeV (the closer to mz the 
better) while m# is driven radiatively to small rather than large values. The TeV-scale top 
squark masses are highly mixed which uplifts rrih to ~ 125 GeV whilst suppressing radiative 
corrections to the scalar potential which influence the values of mh,z- While one expects a 
higgsino-like LSP under conditions of gaugino mass universality, models with non-universal 
gaugino masses allow for a bino-like or wino-like LSP without sacrificing naturalness [33]. 
Mixed axion-higgsino dark matter has been previously calculated in Ref’s [43, 44, 34] while 
the mainly bino-like LSP case is largely excluded due to overproduction of WIMPs [43, 34]. 
Here, we consider the wino-like LSP case which typically yields a thermally-produced wino 
abundance of ~ 0.001 for winos with m^ ~ 100 — 200 GeV (at least an order-of- 

magnitude lower than expectations for a similarly massive higgsino LSP). 

2.8 Two benchmark points 

In order to compute the mixed axion/wino dark matter relic abundance in the SUSY axion 
models, we must specify both the PQ and the MSSM parameters. On the MSSM side, we 
adopt two SUSY benchmark models for illustration. 

The first has been listed as benchmark CSB since it occurs in the rather simple and 
elegant charged SUSY breaking model of Ref. [9]. It is rather similar to the Kallosh- 
Linde [14] benchmark from the study of Ref. [45]. We take the CSB benchmark to be 
illustrative of the large class of models with multi-TeV scalars but with sub-TeV gauginos 
with a wino as LSP. The CSB benchmark model is listed in Table 1. We generate the 
SUSY model spectra with Isajet 7.83 [46]. 
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CSB 

RNSw 

m 0 

72000 

5000 

Mi 

1320 

700 

m 2 

200 

175 

m 3 

-600 

700 

Ao 

0 

-8000 

tan/3 

10 

10 


3000 

200 

rn A 

72000 

1000 

m h 

126.0 

124.3 

m-g 

1924 

1810 

ma, L 

71830 

5101 

m u 

47760 

1478 

m z 2 

635.9 

211.8 

rriy 
&i 

203.2 

114.2 

A-ew 

22830 

10.78 

n s ± d h 2 

0.0020 

0.0015 

Z 1 __ 

a sl (Zip) pb 

6.2 x 10“ 12 

4.3 x 10“ 8 

a SD (Zip) pb 

1.4 x 10" 8 

9.0 x 10~ 4 

{av) „ =0 cm 3 /s 

1.7 x 10“ 24 

1.7 x 10" 24 


Table 1: Masses and parameters in GeV units for two benchmark points computed with Isajet 
7.83 and using m t = 173.2 GeV. 

Along with the CSB benchmark, we adopt a natural SUSY benchmark with a wino 
as LSP. It is taken from Ref. [33] and denoted as RNSw (radiatively-driven natural SUSY 
with a wino LSP). 

3. Brief review of coupled Boltzmann calculation 

To accurately estimate the mixed axion/neutralino dark matter production rate in the 
early universe, it is necessary to evaluate the coupled Boltzmann equations which track 
dark matter number densities and energy densities in an intertwined manner. The exact 
equations used are presented in Ref. [34] and will not be repeated here. In our calculations, 
we use a combination of IsaReD [47] and micrOMEGAs [48] for the evaluation of the wino 
annihilation cross-section ((av)). 

The relevant equations track the following number and energy densities: 

1. neutralino densities including thermal production and production via decays of heav¬ 
ier partices ( e.g. axinos, saxions and gravitinos) followed by possible subsequent 
re- annihilation, 

2. thermally-produced axinos along with axino production via heavy particle decays and 
diminution of axinos due to their decays, 
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3. thermally produced saxions along with diminution via their decays, 

4. thermally-produced gravitinos [49] along with gravitino decay [50], 

5. thermally-produced axions along with axion production via heavy particle decays, 

6. axion production via coherent oscillations (CO) and 

7. saxion production via CO along with saxion decays. 

8. Along with these, we track the radiation density of SM particles. 

For thermal saxion and axion production, it is reasonable to expect annihilation/production 
rates which are similar to axinos. 

The above eight components result in 16 coupled Boltzmann equations: one for the 
number density and one for the energy density of each component. Together with the 
Friedmann equation H = \Jpr/ 3Mj, (where pp is the energy density summed over all 
contributions and Mp is the reduced Planck scale) the Boltzmann equations form a closed 
system which may be solved numerically. 

For the SUSY KSVZ model, the various axino (a — > gg , ZZi and yZ*) and saxion 
branching fractions (s —> gg , gg) can be found in Ref’s [51, 52]. In addition, the model- 
dependent decays s —> aa, aa are effectively parameterized [59] by £ = JT Qi v i / v pq where 
qi are the charge assignments of PQ multiplets and Vi are their vevs after PQ symmetry 
breaking and vpq = y v iQi- We will take ( = 0 or 1 which effectively turns off or 
on saxion decays to axinos/axions [43]. The decay s —> aa augments the LSP abundance 
whilst the decay s —>• aa leads to dark radiation parameterized by the effective number of 
extra neutrinos present in the early universe AV e ff. The Planck Collaboration reported 
N e ff = 3.52^q 45 by the combined data (95%; Planck+WP+highL+Tfo+BAO) [53]. 2 We 
require the upper bound AA r e ff < 1 as a reference value lest too much dark radiation is 
produced. Excluded points with AN e s > 1 are color-coded in our results. 

For the SUSY DFSZ model, axino and saxion decay rates are very different from the 
KSVZ case. While in the KSVZ model axino and saxion decay primarily to gauge bosons 
and gauginos, in SUSY DFSZ then typically a —> Z^ (where cf> = h,H, A), ZiZ, WjW and 
WfH* 1 and s -» pairs of Higgs bosons, vector bosons and electroweak-ino pairs. Complete 
formulae for the DFSZ decay rates are found in Ref. [44], 

The thermal production rates for SUSY KSVZ (which are proportional to Tp) are 
found in Ref’s [55] while thermal production rates for SUSY DFSZ (which are mostly 
independent of Tp) are obtained from Ref’s [56]. We include production of particles via 
both decays and inverse decays [34]: the latter effects are important in SUSY DFSZ where 
saxions and axinos are maximally produced at T ~ m(particle) which leads to a freeze-in 
effect [57] which manifests itself essentially as delayed saxion/axino decays. 

An example of the evolution of various energy densities pi vs. the cosmic scale factor 
R/Ro is shown in Fig. 1 for the SUSY DFSZ model. Ro is taken to be the scale factor at 

2 As this paper was being finalized, this value was updated [54] to N e ff = 3.15 ± 0.23. 
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Figure 1: Evolution of various energy densities vs. scale factor R/Rq for the CSB benchmark case 
in SUSY DFSZ with £ = 1 and other parameters as indicated in the figure. 


the end of reheating (T = Tr). In the figure, f a = 5 x 10 14 GeV while m s = m 3 / 2 = 72 
TeV for the CSB benchmark point. We also take = 40 TeV and £ = 1 so that saxion 
decay to axions is turned on. At R/Ro = 1, the universe is indeed radiation dominated 
(gray curve) while including a thermal population of WIMPs, saxions, axions, axinos and 
gravitinos. It also includes a CO-component of saxions. As R/Ro increases (decreasing 
temperature as denoted by the green dashed line), the oscillating saxion field begins to 
decay- mainly via s —» aa~ so that the population of thermal/decay-produced axions (red 
curve) increases beyond its otherwise thermal trajectory around and below R/Ro ~ 10 4 . 
The neutralino abundance (dark blue) begins to freeze-out around R/Ro ~ 10 5 , but then is 
augmented by decaying CO saxions and also by axinos (which decay slightly after saxions). 
Decaying gravitinos add, but only marginally, to the neutralino abundance around R/Ro ~ 
10 10 . At R/Ro ~ 10 8 , the axion mass turns on and the axion field begins to oscillate 
as non-relativistic matter (brown curve). Also, at R/Ro ~ 10 9 , the neutralinos become 
non-relativistic. Together, the combined neutralino-axion CDM ultimately dominates the 
universe at around R/Ro ~ 10 16 . The ultimate dark matter density is composed of ~ 25% 
wino-like WIMPs and ~ 75% cold axions with a modest-but-not-yet-excluded contribution 
of relativistic axions (AN e ff = 0.68) as dark radiation. 

4. Mixed axion-wino dark matter 

In the following subsections, we compute the neutralino and axion relic abundances for 
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the two benchmark points through numerical integration of the Boltzmann equations as 
discussed in Sec. 3. To gain more general results, we will scan over the PQ scale f a and 
the axino mass which we take to be bounded by 1713 / 2 - 

10 9 GeV < f a < 10 16 GeV, 

0.4 TeV < m~ a < m 3 / 2 , (4.1) 

with m s fixed as m s = m 3 / 2 . In many supergravity models, saxion mass is generated 
by the same operators as those for the MSSM scalars while axino mass is highly model 
dependent and can be much smaller than m 3 / 2 [58, 59]. For this reason, we consider the 
above parameter range for our general analyses. 

For simplicity, we will fix the initial saxion field strength, which sets the amplitude of 
coherent saxion oscillations, to s, = f a (8 S = Si/f a = 1). In addition- for points which 
are DM-allowed (11^ h 2 < 0.12) and obey BBN and dark radiation constraints- the initial 
axion mis-alignment angle is set to the required value such that Q ^ ti 2 + Q a h 2 = 0.12. 
In the SUSY DFSZ case, unlike the SUSY KSVZ model, the bulk of our results do not 
depend strongly on the re-heat temperature (Tr) since the axion, axino and saxion TP rates 
are largely independent of this quantity. Nonetheless, the gravitino thermal abundance is 
proportional to Tr and since gravitinos are long-lived they may affect BBN or WIMP 
abundance constraints if Tr is sufficiently large. In order to avoid the BBN constraints 
on gravitinos, we choose Tr = 10' GeV, which results in a sufficiently small (would-be) 
gravitino abundance. As a result, gravitinos typically do not contribute significantly to the 
neutralino abundance, as discussed above. 

For each of the CSB and RNSw benchmark points, we consider two different cases: 
£ = 0 (saxion decay to axions/axinos turned off) and £ = 1 (saxion decay to axions/axinos 
turned on). We adopt a KSVZ model with SU(2 )r singlet heavy quark states so that the 
axion superfield only has interactions with SU(3) C and U(l)y gauge superfields. We discuss 
the case of SU(2 )r doublet heavy quark states in Sec. 6 for completeness. 

4.1 CSB benchmark in SUSY KSVZ 

4.1.1 £ = 0 case 

In this section, we will examine the CSB benchmark in the SUSY KSVZ case. We start 
with the case where saxion decays into axinos and axions are turned off ( £ = 0). Results 
for this benchmark are shown in Fig. 2, where we plot D^/ 1 2 (blue points) and D a h 2 (green 
points) vs. f a for the scan over parameters defined in Eq. (4.1). In the figure, red points 
violate BBN bounds on late-decaying neutral relics [60] while otherwise the points are BBN 
safe. We also show the measured abundance of CDM by the solid horizontal line. Points 
above this line are excluded by overproduction of dark matter while points below the line 
are allowed. The dashed horizontal grey line denotes the 50% CDM abundance so that 
blue points above this line have WIMP-dominated CDM while green points above this line 
have axion-dominated CDM. 

In Fig. 2, one can see that there are three branches of the neutralino CDM density 
for fa 1$ 10 15 GeV. These branches reflect three regions of axino mass. The uppermost 
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Figure 2: The wino-like WIMP (blue) and axion (green) relic densities from a scan over SUSY 
KSVZ parameter space for the CSB benchmark case with £ = 0. The grey dashed line shows the 
points where DM consists of 50% axions and 50% neutralinos. 


branch corresponds to the case of < rn^. In this case, axinos decay only into Z\ plus 
SM particles. Since the axion sector does not have a direct coupling to a SU(2 )l gauge 
supermultiplet, axino decays into Z\ (mostly wino-like) happen only through the bino-wino 
mixing, which is very tiny in the MSSM. Therefore, for this branch the axino decay occurs 
well-after the neutralino freeze-out, enhancing the neutralino abundance well above the 
measured CDM density for all values of f a . Moreover- for f a > 10 10 GeV- BBN constrains 
the model due to the long-lived axino. 

The middle branch corresponds to < rrig. In this region axinos can decay 

directly into Z 2 . Since Z 2 is mostly bino-like and axinos directly couple to B through the 
U(l)y anomalous coupling, their life-time is much shorter than in the rria < m z 2 case. 
Although the axinos decay after neutralino freeze-out for all f a , the neutralino density still 
is smaller than the observed CDM density for f a < 5 x 10 10 GeV. Hence, both axion- 
dominated or neutralino-dominated dark matter scenarios are possible in this region. For 
fa > 5 x 10 12 GeV, all points in the rria < rrig branch are excluded by BBN. 

The lowermost branch corresponds to > rrig. In this region, axinos can decay to 
gluinos through the SU(3) C anomaly coupling so that the axino life-time becomes much 
shorter than the previous two cases. For f a < 10 12 GeV, axinos decay before neutralino 
freeze-out in the bulk of this parameter region, so the neutralino CDM density takes its 
standard thermal value ~ 0.002. In the case where the axino mass is close to the gluino 
mass, however, axinos can decay after neutralino freeze-out and augment the WIMP abun¬ 
dance. As f a increases, axinos more often decay after freeze-out and hence increasingly 
augment the neutralino relic density. By f a ~ 2 x 10 12 GeV, axinos always decay af¬ 
ter freeze-out and always augment the neutralino abundance. Despite the enhancement 
of the neutralino abundance, there are points where the DM is axion-dominated up to 


10 









f a (GeV) 


Figure 3: The wino-like WIMP and axion relic densities from a scan over SUSY KSVZ parameter 
space for the CSB benchmark case with £ = 1. The grey dashed line shows the points where DM 
consists of 50% axions and 50% neutralinos. 


/„ ~ 6 x 10 13 GeV. 

For fa > 10 15 GeV, the contribution to the WIMP abundance is mostly from CO- 
produced saxion decays; these augment the abundance for larger f a since the saxion CO 
production rate increases with f a . On the other hand, the contribution from TP axinos is 
highly suppressed for large f a , since the axino thermal production decreases with f a . Once 
the TP axino abundance becomes negligible (f a > 5 x 10 14 GeV), the LSP relic abundance 
becomes independent of the axino mass and all the branches discussed above collapse into 
a single line, as seen in Fig. 2. 

4.1.2 £ = 1 case 

In Fig. 3, we show 0^ Jh 2 vs. f a for the same CSB benchmark point but now where 
saxion decays into axinos and axions are allowed: £ = 1. For the lower (f a < 10 14 GeV) 
range, saxion decays have a smaller impact on the neutralino abundance and the results are 
similar to the CSB/KSVZ £ = 0 case. For higher f a values, CO-produced saxions become 
important and since s —> aa and aa decays are now allowed, there is a large injection 
of relativistic axions. For f a > 4 x 10 14 GeV, we see brown points which produce too 
much dark radiation— AW e ff > 1“ and are excluded. There is also a broad band of blue 
(BBN-allowed) and red (BBN-excluded) points at large f a ~ 10 15 GeV with very high 
n z 1 h ‘ 2 ~ 1 — 100 where the additional neutralino abundance arises from s —» aa decays. 
The lower disjoint narrow band at f a > 10 14 GeV and G^/j 2 ~ 0.1 — 1 occurs for points 
where > m s / 2, so s —> aa is kinematically forbidden. 

Finally we point out that, unlike the £ = 0 case, the extremely large f a region 
( f a > 10 15 GeV) still shows a dependence on the axino mass: this is responsible for the 
distinct branches. Although thermal production of axinos is neglegible in this regime, axi- 
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Figure 4: The wino-like WIMP and axion relic densities from a scan over SUSY DFSZ parameter 
space for the CSB benchmark case with £ = 0. The grey dashed line shows the points where DM 
consists of 50% axions and 50% neutralinos. 


nos are non-thermally produced from saxion decays and can influence the final neutralino 
abundance. 

4.2 CSB benchmark in SUSY DFSZ 
4.2.1 £ = 0 case 

In this section, we will examine the CSB benchmark in the SUSY DFSZ case. As before, 
we start with the £ = 0 case, shown in Fig. 4. The first noteworthy point is that the 
large /r value enhances the saxion (axino) decay rate to Higgs (higgsinos). As a result the 
saxion and axino lifetimes are suppressed and the entire f a range is BBN safe. Unlike the 
KSVZ case, there are two branches for neutralino CDM density, since in the DFSZ case, 
the axino decay is determined by the ^-term interaction. The upper branch corresponds to 
rria < m z 3 ~ M with higgsino-like Z 3 . The axino decay into Z\ or Z 2 can be through wino- 
higgsino or bino-higgsino mixing, so it is normally suppressed by (mz/n) 2 - For f a < 3x 10 10 
GeV, axinos decay before neutralino freeze-out, and thus the neutralino density takes its 
standard value. For / a > 3 x 10 10 GeV, axinos tends to decay after neutralino freeze- 
out so the neutralino density gradually increases as f a increases. In most of parameter 
space, axions constitute the bulk of dark matter, but wino-like neutralinos can be the 
dominant dark matter in the region of 10 12 GeV< f a < 10 13 GeV. By f a > 10 13 — 10 14 
GeV, the neutralino density is typically larger than the measured CDM result so the 
parameter/model choices would be excluded. 

The lower branch corresponds to rria > m z 3 - Due to its large interaction, the axino 
tends to decay before neutralino freeze-out for / a < 3 x 10 12 GeV. Therefore, the neutralino 
relic abundance is usually fixed at its thermally-produced value for much of the lower range 
of f a . Once f a > 10 13 GeV, the neutralino abundance is always enhanced due to decays 
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Figure 5: The wino-like WIMP and axion relic densities from a scan over SUSY DFSZ parameter 
space for the CSB benchmark case with £ = 1. The grey dashed line shows the points where DM 
consists of 50% axions and 50% neutralinos. 


of axinos and saxions. Still, the CDM abundance tends to be axion-dominated for f a < 
2 x 10 14 GeV. For higher f a there is a short interval where wino-like WIMPs can dominate 
the DM abundance. Finally, for f a > 6 x 10 14 GeV, WIMP CDM is always overproduced. 
We also point out that for very large f a values, as in the KSVZ £ = 0 scenario, the 
thermal production of axinos is neglegible hence the neutralino relic abundance becomes 
independent of f a . 

4.2.2 £ = 1 case 

For the CSB benchmark with SUSY DFSZ and £ = 1, the results are shown in Fig. 5. The 
low f a behavior of the plot is similar to the CSB/DFSZ case with £ = 0: the CDM density is 
dominated by axions. For higher f a values, where CO-produced saxions become important, 
the saxion lifetime is shortened by the additional contributions from s — > aa, aa decays. 
However, most of the points for /„ > 5 x 10 14 GeV are forbidden due to overproduction 
of dark radiation. The lower blue-brown band at f a 10 14 — 10 15 GeV occurs when 
rria > rn s /2 so that additional WIMP production from s —> aa is dis-allowed. 

4.3 RNSw benchmark in SUSY KSVZ 

In this subsection we examine dark matter production in the SUSY RNSw benchmark case. 
The RNSw benchmark model has values m s = m o = m 3 / 2 = 5 TeV which is far smaller 
than that of the CSB benchmark so that saxions (and also axinos since we take their mass 
to be bounded by m 3 / 2 ) are typically much longer-lived than in the CSB case. 

4.3.1 £ = 0 case 

In Fig. 6, we plot Qg a h 2 vs. f a for the SUSY KSVZ case using the RNSw benchmark point 
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Figure 6: The wino-like WIMP and axion relic densities from a scan over SUSY KSVZ parameter 
space for the RNSw benchmark case with £ = 0. The grey dashed line shows the points where DM 
consists of 50% axions and 50% neutralinos. 


with £ = 0. In this case, rria is always larger than to ^, so there are only two branches for 
the neutralino density: rria < rrig and rria > rrig. Long-lived axinos are already augmenting 
the neutralino relic density even at f a values as low as 10 9 GeV. As we move to higher 
f a values, the axinos and saxions are longer-lived, thus contributing even more to the 
WIMP abundance. For f a values > 5 x 10 12 GeV, the model is already excluded due to 
overproduction of WIMPs. The fl^/i 2 points reach even higher values as f a increases until 
/ a ~ 3 x 10 13 GeV. For 3 x 10 13 GeV< f a < 2 x 10 14 GeV, the axino contribution decreases 
due to suppression of the thermal production. For / a > 2 x 10 14 GeV, then CO-produced 
saxions decay into gluino pairs and tend to augment the WIMP abundance. However, this 
is inconsequential since the model already overproduces WIMP dark matter. A large BBN 
forbidden region occurs, but it is already in the WIMP-overproduction region so adds no 
further constraints. 

4.3.2 £ = 1 case 

For the RNSw benchmark case in SUSY KSVZ with £ = 1, as shown in Fig. 7, the low f a 
behavior of is very similar to the £ = 0 case, since at low f a saxion production is not 

very relevant and saxions decay well-before neutralino freeze-out. For / a >4x 10 12 GeV, 
as in the £ = 0 case, the model over-produces WIMPs and is excluded. At even larger 
values of f a , in the DM-excluded region, the £ = 1 case begins to differ from £ = 0. An 
additional branch of flj^/i 2 appears: the lowermost branch swings downward due to the 
suppressed axino and saxion TP as in the £ = 0 case, but never reaches the DM-allowed 
line. It is nonetheless also excluded by overproduction of dark radiation. The upper (also 
excluded) two branches occur where m s > 2 rria so that CO-production of saxions keeps 
increasing the WIMP abundance. This region is also excluded by the BBN constraint from 
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Figure 7: The wino-like WIMP and axion relic densities from a scan over SUSY KSVZ parameter 
space for the RNSw benchmark case with £ = 1. The grey dashed line shows the points where DM 
consists of 50% axions and 50% neutralinos. 


late-decaying saxions followed by axino cascade decays. 

4.4 RNSw benchmark in SUSY DFSZ 

In this subsection, we examine the RNSw benchmark model in the SUSY DFSZ model. 
Since /r(RNSw) <C ^(CSB), saxions in the RNSw/DFSZ case will tend to be longer lived. 

4.4.1 £ = 0 case 

In Fig. 8 we show RNSw in the SUSY DFSZ case with £ = 0. In this case, rria is always 
larger than fi, so there is no region corresponding to the upper branch in Fig. 4. For low / a , 
in contrast to RNSw in the SUSY KSVZ case, the axino lifetime is smaller and the WIMP 
abundance remains at its thermally-produced value for /„ < 5 x 10 10 GeV. For higher f a 
values, the WIMP abundance is augmented by axino and saxion decays after freeze-out. 
Ultimately, the model over-produces WIMPs for f a > 10 13 GeV. The model tends to be 
axion-dominated for f a < 6 x 10 12 GeV and WIMP dominated for a narrow range of f a just 
beyond this value until WIMP overproduction is reached and the model becomes excluded. 
This is in contrast to the CSB benchmark with DFSZ and £ = 0, where the allowed region 
extends to f a ~ 5 x 10 14 GeV since saxions and axinos are shorter-lived due to much larger 
masses and stronger interactions (/r(RNSw) <C /r(CSB)). 

4.4.2 £ = 1 case 

In Fig. 9 we show results for the RNSw benchmark in SUSY DFSZ with £ = 1. While 
the low f a behavior is similar to the results from the £ = 0 case, the high f a behavior is 
different. The decays s —X aa and s —x aa allow the saxion to decay more quickly than in 
the £ = 0 case for a common value of f a ■ Thus, the DM-allowed region extends to larger f a 
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Figure 8: The wino-like WIMP and axion relic densities from a scan over SUSY DFSZ parameter 
space for the RNSw benchmark case with £ = 0. The grey dashed line shows the points where DM 
consists of 50% axions and 50% neutralinos. 
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Figure 9: The wino-like WIMP and axion relic densities from a scan over SUSY DFSZ parameter 
space for the RNSw benchmark case with £ = 1. The grey dashed line shows the points where DM 
consists of 50% axions and 50% neutralinos. 


values: in this case up to f a ~ 10 14 GeV. For these high f a values, the relic density band 
again splits into two branches: one with heavy axinos (lower-branch), where s —> aa is 
closed, and one with light axinos (upper branch), where s —>• aa is open, thus augmenting 
the WIMP abundance. The points with / a > 2 x 10 14 GeV tend to be doubly-excluded by 
overproduction of WIMPs and by overproduction of dark radiation. 
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Figure 10: Plot of thermally-produced neutralino abundance vs. m 3 / 2 along the CSB 

model line with tan/3 = 10 and /r = 3 TeV. 


5. Dependence of mixed axion-wino abundance on sparticle mass spectra 

In the previous sections we have investigated the DM-allowed range of f a for two SUSY 
benchmark models with wino-like LSPs, f a and as free parameters and m s = mo- In this 
section, we investigate how our results might change as a function of the MSSM spectrum. 
To explore this issue, we extend our two benchmark points into model lines in the MSSM 
sector. For brevity, we consider here only the DFSZ model which provides a solution 
to the SUSY fj, problem- to see the impacts of axino/saxion production/decays on the 
CDM density. Actually, even in the presence of late-decaying axinos and saxions, the most 
important factor that determines the WIMP abundance is the WIMP-WIMP annihilation 
cross section since the augmented density is determined mainly by annihilation cross section 
evaluated at the heavy particle decay temperature: this is the case of so-called WIMP re¬ 
annihilation after non-thermal WIMP production from heavy particle decay [51]. For this 
reason, the behavior of our plots is similar for both DFSZ and KSVZ models, and so we 
will show only the DFSZ case and then briefly comment on the KSVZ case. 

5.1 CSB model line 

For the CSB benchmark, we will now allow 7713/2 to vary while keeping tan/3 fixed at 10 
with [i = 3 TeV and m,A = 7773 / 2 . For the CSB model-line, we require m 3 / 2 > 32 TeV 
so the mass of the lightest wino-like chargino is always above the limit > 91.9 GeV 
established from LEP2 searches. The upper limit on m 3 / 2 occurs at ~ 115 TeV where the 
predicted value of rrih climbs above 128 GeV. Here, we allow for an expected theory error 
in the Isasugra calculation of rrih at about ±2.5 GeV. 

We show the thermally-produced neutralino abundance for the CSB model line in 
Fig. 10. Here, we see that D~ p /i 2 ranges from around 0.0007 at the lower limit to about 
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Figure 11: Plot of thermally- and non-thermally-produced neutralino abundance vs. f a 

along the CSB model line in the DFSZ £ = 1 case for a light (m 3 / 2 = 32 TeV, blue envelope) and 
heavy (m 3 / 2 = 115 TeV, green envelope) CSB mass spectrum where rria ranges from 400 GeV up 
to m 3 / 2 . 


0.005 at the upper limit as compared to 0.002 for the CSB benchmark. Roughly speaking, 
the thermally-produced wino abundance will provide either more or less room in general 
for non-thermally produced winos and axions. 

In Fig. 11, we show the value of which is produced from the coupled Boltzmann 

calculation of mixed axion-wino CDM versus f a for the minimal and maximal values of 
777 , 3/2 which are allowed along the CSB model line. The blue curves provide the calculated 
envelope of values for the lower limit of 7773/2 I'N-' 32 TeV. At low f a , 12^/i 2 lies at the TP- 
value ~ 0.0006 since thermally-produced axinos always decay before neutralino freeze-out. 
As f a climbs above ~ 10 11 GeV, then the lighter axinos start decaying after neutralino 
freeze-out whilst the heavier axinos still decay before freeze-out. The region between the 
two blue curves shows the range of H^/t 2 which is generated for 0.4 TeV < rria < 32 
TeV. We see that values of f a up to ~ 10 15 GeV are dark-matter-allowed for very heavy 
axinos. However, at values of f a > 5 x 10 14 GeV, then too much dark radiation is produced 
from s —^ aa decays in addition to WIMP overproduction so that the parameter space is 
doubly-excluded. 

The heavy end of the CSB model line 7773/2 = 115 TeV is shown by the envelope of 
green curves. For the light axino with rria = 421 GeV, the thermally-produced value of 
12 2 h 2 ~ 0.003 is obtained only for the short range of f a < 3 x 10 9 GeV. For this heavy 
CSB spectra, the gauginos are all sufficiently heavy and thus the axino can decay only into 
Z\ so that the axino lifetimes are much longer than that in the case for light spectra. The 
upper range of the green envelope comes from light axino masses where rria = 421 GeV 
is the threshold for a ZZ\ decay hence augmenting neutralino density at low f a , while 
the lower envelope is established by the heaviest axino mass values. For the upper part 
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Figure 12: Plot of thermally-produced neutralino abundance vs. M 2 along the RNSw 

model line. 


of the envelope, the red points denote the on-set of BBN bounds on late decaying saxions 
as ruling out / a > 5 x 10 13 GeV. For the lower part of the envelope, with axino masses 
ranging to 32 TeV, then values of f a up to 5 x 10 14 GeV are possible. 

In the case of KSVZ model, the spectrum dependence is similar to the DFSZ model. 
For the light spectum ( 777 , 3/2 = 32 TeV), the neutralino abundance tends to be smaller due 
to its large annihilation cross section. For the heavy spectrum ( 7773/2 = 115 TeV), the cross 
section becomes larger, so the neutralino abundance becomes smaller. Nevertheless, the 
allowed range of f a for the heavy spectrum is slightly larger than that for light spectrum 
since the saxion mass is larger (m s = 777,3/2) so that its decay can occur earlier. 

5.2 RNSw model line 

For the RNSw benchmark, we will instead allow the GUT scale SU(2)^ gaugino mass M 2 
to vary while keeping iuq = 5 TeV, 7773/2 = 700 GeV, Aq = —8 TeV and tan/3 fixed 
at 10 with /x = 200 GeV and 777,4 = 1 TeV. For the RNSw model-line, the lower limit 
on M 2 is again set by the limit from LEP2 searches for wino-like charginos. The upper 
limit on M 2 < 250 GeV is set from simply requiring a wino-like LSP: for higher M 2 
values, the lightest neutralino becomes increasingly higgsino-like, a case which was shown 
in Ref’s [44, 34], The naturalness value A ew remains fixed at around 10 since varying M 2 
hardly affects it [33]. 

We show the thermally-produced neutralino abundance for the RNSw model line in 
Fig. 12. The lower range of H~ p /i 2 occurs at 0.0012 for M 2 140 GeV. The maximal value 
reaches up to ~ 0.004 before entering the higgsino-like LSP region. Since the thermally- 
produced wino abundance increases with M 2 , the allowed enhancement from non-thermal 
production decreases as M 2 increases. Furthermore, since the non-thermal production 
(from saxion and axino decays) grows with f a , we expect the maximum allowed value for 
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Figure 13: Plot of upper limit of f a allowed from the RNSw model line in the DFSZ case with 
£ = 1 versus M 2 for mj = 400 GeV and ma = 5 TeV. 


f a to decrease as M 2 increases. This is shown in Fig. 13, where we plot the upper limit 
on f a , denoted as /*, versus M 2 along the RNSw model line in the DFSZ £ = 1 case. The 
axino mass is rria = 0.4 TeV (green dots) or = 5 TeV (black dots). The upper limit 
comes only from the overproduction of WIMPs in this case since violation of BBN bounds 
and overproduction of dark matter occurs at higher f a in this model. 

In the KSVZ model, on the other hand, axinos are longer-lived than in the DFSZ 
model, so the allowed range of f a is smaller than that in the DFSZ case. As we have 
seen in Fig. 7, only a small region, f a < 0(1O 10 ) GeV, is allowed for < rrig, while 
f a < 0(1O 12 ) GeV is allowed for = 5 TeV. 

6. Conclusion 

In this paper we have examined mixed axion/wino cold dark matter production in two 
SUSY benchmark models with a wino as LSP. The first- labeled as CSB- is typical of a 
variety of models (PeV-SUSY, some split SUSY variations, KL, PGM, spread SUSY) with 
a thermally-underproduced wino-like WIMP abundance. The second, labeled as RNSw, is 
a model with radiatvely-driven naturalness but with a wino-like rather than a higgsino-like 
LSP. Our calculation of mixed axion/wino dark matter production stands in contrast to 
the more commonly examined case of non-thermal WIMP production due to late decaying 
moduli fields [4, 19, 20]. We find it a more appealing method for augmenting the dark 
matter abundance since it also provides a solution to the strong CP problem and- in the 
case of SUSY DFSZ- provides for a solution to the SUSY // problem. 

We have presented results for the wino-like WIMP abundance and axion abundance 
as a function of the axion decay constant f a and the axino mass rria■ In the bulk of 
the parameter space, WIMPs are thermally under-produced at low and intermediate f a 
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Figure 14: Range of f a which is allowed in each PQMSSM scenario for the CSB and RNSw 
benchmark models. Darker-shaded regions indicate the range of f a where 9i > 3 which might be 
considered unnatural. We also show the f a range which is expected to be probed by the ADMX 
experiment in the next few years. 


values (~ 10 9 — 10 11 GeV) so that the DM abundance tends to be axion-dominated. This 
has important consequences for direct and indirect WIMP detection experiments since it 
anticipates a greatly reduced local abundance of WIMPs and hence diminished prospects 
for wino-like WIMP detection. This can actually allow for wino-like WIMP dark matter 
to evade the recent Fermi [61] searches for gamma ray emission from dwarf-spheroidal 
galaxies since in this case the expected event rate is expected to be reduced by a factor 

(n w h 2 / o.i2) 2 . 

A grand overview of our results is presented in Fig. 14, where we show the allowed 
range of f a as a bar for each of the two benchmark points and for each of the eight SUSY 
PQ models considered. For all the models, no GUT scale values of f a (f a ~ 10 16 GeV) are 
allowed. This is due to the rather large value of m s ~ 7713/2 in our benchmark models. In 
these cases, saxions always decay to SUSY particles and no entropy dilution of WIMPS 
and axions is possible (see Refs. [44, 34] for more details). 

In addition to the SUSY KSVZ case with SU(2)^ singlet heavy quark states which 
has been presented here, we have also investigated SUSY KSVZ models including SU(2)i 
doublet heavy quark states so that the axion superfield has couplings with SU(2)^ gauge 
superfields. In the case of doublet heavy quarks, axino decays to the wino-like neutralino 
are not suppressed, even for rria < m z 2 ■ Therefore, there is no separate branch like the 
uppermost one in Fig. 2 and 3 and thus there are only two branches determined by m,g. 
The basic features of plots with doublet heavy quarks are similar to the case with singlet 
heavy quarks since the dominant axino decay mode is into gluinos for both cases. The 
allowed range of f a values is extended only slightly for doublet KSVZ heavy quarks as 
compared to the case of singlet heavy quarks shown in this paper. 
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For sufficiently heavy axinos, all models shown in this paper are DM-allowed for the 
lower range of f a 10 9 - 10 12 GeV, since WIMPs are underproduced. In these cases, the 
remaining abundance is made up of axions. Even though one might expect a low axion 
abundance at low f a in the case where the initial mis-alignment angle is Oi ~ 0(1), due 
to anharmonicity effects the necessary axion abundance can always be obtained by taking 
Qi ~ 7 r [31]. In this case, one might wonder about fine-tuning of the axion abundance such 
that the axion fields sits atop the peak of its potential. Thus, for cases where 9{ > 3, we 
shade these regions as darker in Fig. 14. The non-shaded regions may be more natural as 
far as the expected initial axion field value goes. 

We should note that for KSVZ models, regions with 9i < 3 at low f a occur only if the 
wino LSP constitutes more than ~ 90% of total CDM density since axion CO-production 
is very low for f a < 10 10 . Then, the CSB benchmark in the SUSY KSVZ model most 
naturally allows for the lowest f a values while the CSB benchmark in the DFSZ model 
allows for the highest f a values. The range of f a values obtained for the RNSw benchmark 
is more constrained than the CSB case. The upper bounds on f a for the two benchmark 
models are well-maintained even when the points are extended to model lines, as was shown 
for the DFSZ £ = 1 case in Sec. 5. 

Finally, we denote the range of f a values which are expected to be probed in the 
next few years by the ADMX experiment [62], The values shift between KSVZ and DFSZ 
models since the domain wall number IVdw = 1 for KSVZ and 6 for DFSZ and m a — 
0.62 eV[10' GeV/(/ a /-/VDw)]- We also note that a possible ADMX technique of open 
resonators [63] may allow even lower values of f a to be probed in the future. 
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